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The binding energies of ^®Si, ^^S, and "^^Cl have been measured with a relative uncertainty 
< 0.59 X 10~® using a flat-crystal spectrometer. The unique features of these measurements are 
1) nearly perfect crystals whose lattice spacing is known in meters, 2) a highly precise angle scale 
that is derived from first principles, and 3) a gamma-ray measurement facility that is coupled to 
a high flux reactor with near-core source capability. The binding energy is obtained by measuring 
all gamma-rays in a cascade scheme connecting the capture and ground states. The measurements 
require the extension of precision flat-crystal diffraction techniques to the 5 to 6 MeV energy region, 
a significant precision measurement challenge. The binding energies determined from these gamma- 
ray measurements are consistent with recent highly accurate atomic mass measurements within a 
relative uncertainty of 4.3 x 10~^. The gamma-ray measurement uncertainties are the dominant 
contributors to the uncertainty of this consistency test. The measured gamma-ray energies are in 
agreement with earlier precision gamma-ray measurements. 



I. INTRODUCTION 

Nuclear binding energy measurements are of interest 
because they are accurately related to atomic mass mea- 
surements. This relationship provides a means to check 
the results in one precision measurement field against re- 
lated results obtained in another precision measurement 
field. Because the experimental techniques used in the 
two fields are very different, this check has the potential 
to reveal systematic errors associated with either mea- 
surement. 

The synergism between binding energy and atomic 
mass measurements can be demonstrated by considering 
a typical neutron capture reaction n + '^X ^+^X -I- 7's 
which leads to the following equation involving atomic 
masses and the binding energy. 



m{n) + m('^X) = m{^+^X) + Sn 



(1) 



Atomic masses m are measured in atomic mass units 
while the binding energy of ^"""^X, Sn, is obtained from 
gamma-ray wavelengths measured in meters. The bind- 
ing energy in meters can be converted to atomic mass 
units using the molar Planck constant, NaH, divided by 
the speed of light, c 0. This combination of constants is 
known with a relative uncertainty of 6.7 x 10~^, an ac- 
curacy which does not limit the test implied by Eq. 
given the presently available accuracy in atomic mass and 
binding energy measurements 01 ■ 

New precision measurements of the ^^Si, '^^S, and "^^Cl 
binding energies have been made using a flat crystal spec- 
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trometer. This spectrometer measures the wavelengths 
of the gamma-ray photons using crystals whose lattice 
spacings are known in meters and an angle scale that is 
derived from first principles. Thus, the measured wave- 
lengths are on a scale consistent with optical wavelengths 
and the SI definition of the meter. The binding energies 
of these three nuclei are in the 8.5 to 8.6 MeV range and 
are obtained by measuring lower energy lines that form 
a cascade scheme connecting the capture and ground 
states. For all three nuclei, the cascade scheme with 
the most intense transitions includes a gamma-ray with 
energy > 4.9 MeV. Such high energies present a signifi- 
cant measurement challenge for gamma-ray spectroscopy 
because the Bragg angles are < 0.1° and the diffracted 
intensity is rather small (a few s~^ or less). 



II. THE ^''SI, 3^S, AND ^"^CL DECAY SCHEMES 

In Fig. n we show partial decay schemes for ^^Si, ^'^S, 
and ^^Cl containing the transitions that were measured 
in these binding energy determinations. The values in 
parentheses are the number of gamma-rays emitted per 
100 captures 0]. The reactions, the thermal neutron cap- 
ture cross sections, and the nominal energies of the mea- 
sured gamma-rays are listed in Table|l| 4]. Because ^^Cl 
has the largest thermal neutron capture cross section of 
any light nuclei, the '^^Cl binding energy was measured 
first. The experience gained in increasing the crystal re- 
flectivity and lowering the background during the ^^Cl 
measurement proved to be very valuable in the measure- 
ment of the weaker ^^Si and '^^S transition energies. In 
addition, a larger than expected dependence of the angle 
calibration on the environment was noted during the CI 
measurements which led to more frequent angle calibra- 
tions during the ^^Si and ^"^S measurements. The ^^Si 



2 



33, 



36, 



8472.9 keV \ 8641 .6 keV ^ 8579.8 keV 



3539.1 
(68.0) 



4934.4 
(62.7) 



5421 .0 
(59.1) 



2379.6 
(44.5) 

841.0 (75.6) 



1951.2 
(21.7) 



6111.5 
(20.0) 



517 .1 (18.5) 



1164.9* (19.9) 
786.3* (9.6) 



FIG. 1: Partial decay schemes for ^^Si, ^^S, ^^Cl showing the 
transitions that were measured in this study. The numbers 
in parentheses are the number of gammas per 100 neutron 
captures. The transitions marked with asterisks in CI were 
not used to deduce binding energies. 



measurement is particularly difficult because the 4934 
keV line is emitted while the nucleus is recoiling follow- 
ing the emission of the 3539 keV line. Thus, the 4934 
keV profile is significantly Doppler broadened which de- 
creases the accuracy with which the Bragg angle can be 
determined. In each of these nuclei there are other less 
intense transitions that connect the capture and ground 
states. However, the measurement uncertainty of these 
weaker lines will be so large that their contribution to 
the binding energy determination will be small. As the 
stability of the spectrometer improves and techniques to 
measure weaker lines are developed, these weaker lines 
may be used for future binding energy determinations. 



III. EXPERIMENT 

The measurements were made at the Institut Laue- 
Langevin (ILL) using the GAMS4 flat crystal spectrom- 
eter. A detailed description of this facility is available 
to the interested reader in Ref. 0. The discussion of the 
experiment given here will be limited to those aspects 
that are peculiar to these measurements. This spec- 
trometer is located on the reactor floor at the exit of 
a through-tube that has facilities to transport and hold 
sources next to the reactor core. Five beamtime allo- 
cations were devoted to these binding energy measure- 
ments. Chlorine was measured in February 1997 and 
September 1997. Sulfur was measured in September- 
October 1998 and April-May 1999. Silicon was measured 
in October-November 2000. The first measurement cy- 
cle (February 1997) served as a proof of the capability 
of the facility to determine binding energies. In retro- 
spect, it was necessary to exclude all data taken during 
this cycle from the final results because the experimental 
conditions were not sufhciently controlled. 



The stability of the angle calibration of the spectrom- 
eter has been a particularly troublesome experimental 
problem. Throughout the extended period of the mea- 
surements reported here, the angle calibration of the 
spectrometer has been measured many times, both dur- 
ing and between the binding energy measurements. How- 
ever, at the start of these measurements we were not 
aware of the dependence of the angle calibration on the 
humidity and did not perform calibrations as frequently. 
From 1998 forward, more frequent calibrations have pro- 
vided a means to obtain more accurate angle calibra- 
tions. Following considerable data analyses aimed at 
mining the angle calibrations and gamma-ray measure- 
ments for maximum information, we have reached the 
following conclusions concerning the angle calibration: 
a relative uncertainty due to the angle calibration of 
w 0.4 X 10^^ should be included in each measurement 
period. More details concerning the angle calibration are 
given in Sec. IIVI 

At the start of the experiment, it was our intention to 
measure the 786 keV and 1165 keV lines in CI to provide 
a cascade cross-over verification of the 1951 keV line. 
Because of beamtime restrictions and angle calibration 
difhculties, this cascade cross-over verification was not 
realized. However, a very limited data set for the 786 
keV line was recorded and used to provide a value for 
the wavelength of this transition. 



A. Sources 

The source handling mechanism can accommodate 
three sources that are positioned one behind the other 
on the beam axis, a line connecting the source region 
and the center of the axis of rotation of the first crystal. 
The source material was placed in thin-walled graphite 
holders that are supported on "V"s for precise position- 
ing. Two sizes of graphite holders were used. For the two 
chlorine and first sulfur measurements, the inside volume 
of the source holder was 2 mm x 50 mm x 25 mm with 
the 2 mm x 50 mm surface facing the spectrometer. For 
the second sulfur and silicon measurement, the inside vol- 
ume was 3.5 mm x 50 mm x 25 mm with the 3.5 mm x 
50 mm surface facing the spectrometer. The neutron flux 
at the source position is w 5.5 x 10^'* cm^^s^^. In Ta- 
ble ^1 the sources, the source masses, and the estimated 
activities for each of the measurements are given. The 
sources that are used are extremely active to compensate 
for the small effective solid angle (high resolution) of the 
spectrometer. 



B. Spectrometer 

The critical component of the gamma-ray facility is a 
double flat crystal spectrometer that has three unique ca- 
pabilities that are very important for the accurate mea- 
surement of pico-meter wavelengths. First, the diffrac- 
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TABLE L Reactions, cross-sections, and nominal energies associated with the binding energy determinations. 



Nuchde 


Reaction 


Cross-section 


Nominal 


Energies measured 






(x 10-28 in2) 


binding energy 
(keV) 


(keV) 


29Si 


n + 2«Si^29si + ^ 


0.177 


8473 


3539, 4934 


33s 


n + ^^S -» ^='8 + 7 


0.53 


8642 


841, 2380, 5421 


36(^j 




43.6 


8580 


517, 786, 1951, 6112 



TABLE IL Sources, source masses, and estimated activities for each of the measurements. 



Nuclide 


Measurement 
dates 


Source 


Source mass 
(g) 


Estimated activity 
(x 10'=^ Bq) 


2^Si 


Oct-Nov 2000 


Si single crystal 


13.4 


2.6 


33s 


Sept-Oct 1998 


ZnS single crystal 


8.1 


1.4 


33s 


April-May 1999 


ZnS polycrystal 


16.8 


2.9 


36 Q 


Feb 1997 


NaCl 


4.5 


85.4 


36 Qi 


Sept 1997 


NaCl 


4.5 


85.4 



tion crystals are highly perfect specimens whose lattice 
spacing is measured on a scale consistent with the SI def- 
inition of the meter. Second, the diffraction angles are 
measured with sensitive Michelson angle interferometers 
which are calibrated using an optical polygon The 
angle calibration is based on the fact that the sum of the 
external angles of the polygon equals 360°. Third, the 
gamma-ray beam collimation is sufficient to permit the 
measurement of very small diffraction angles (< 0.05°). 
The GAMS4 facility is a precision metrology laboratory 
with the usual attention to vibration isolation, tempera- 
ture control, and environmental monitoring. 



Because descriptions of the profile recording and the 
angle calibration that follow assume some knowledge of 
the angle measuring system of the spectrometer, a few de- 
tails concerning the angle interferometers are provided. 
Each of the two Michelson angle interferometers con- 
tains two corner cube retro-reflectors that are rigidly at- 
tached to the crystal rotation table. As the crystal ro- 
tates, the path length of one arm of the interferometer in- 
creases while the path length of the other arm decreases. 
The angles are measured in whole and fractional inter- 
ferometer fringes with 1 fringe ~ 7.8 x lO^^rad. An 
interferometer fringe can be divided into ~ 1000 parts 
7.8 X IQ-^^rad). The arm that supports the retro- 
reflectors is made out of low expansion invar in order to 
reduce the temperature dependence of the angle interfer- 
ometer. Nevertheless, the temperature of the corner cube 
arm must be taken into account when converting inter- 
ferometer fringe values into angles (see Sec. lIV|) . Because 
the interferometers are in the laboratory environment, all 
angle fringe measurements must be corrected to standard 
pressure, temperature, and humidity conditions. 



C. Crystals 

The gamma-rays measured in this study were 
diffracted by nearly perfect silicon crystals used in trans- 
mission geometry. All of the crystals were cut so that 
the (220) family of planes was available for diffraction and 
have a shape and mounting identical to the crystal shown 
in Fig. 7, Ref.0. Because of the large spread in energy of 
the gamma-rays, two different sets of crystals were used. 
The first set, called ILL2.5, consisted of two crystals of 
nearly equal thickness (« 2.5 mm), and are the same 
crystals used for the measurement of the deuteron bind- 
ing energy Q. The raw material for these crystals was 
obtained from the Solar Energy Research Institute 0]- 
The second set, called ILL4.4&6.9, consisted of two crys- 
tals of unequal thickness (~ 4.41 mm and w 6.95 mm) 
manufactured from raw material obtained from Wacker 
f^. The ILL2.5 crystals were used for the lower energy 
lines (S 841 keV and 2380 keV and the CI 517 keV, 786 
keV, and 1951 keV lines) and initial measurements of the 
CI 6111 keV line. The ILL4.4&6.9 crystals were used for 
the higher energy lines (Si 3539 keV and 4934 keV, S 
5421 keV, and CI 6111 keV lines). The dependence of 
integrated reflectivity on crystal thickness and energy is 
discussed in more detail in Section 11, Ref. 0. Each line 
was measured in at least two orders that were chosen on 
the basis of high reflectivity and small diffraction width. 
In Table imi the crystal configurations that were used for 
each energy are given along with the nominal Bragg an- 
gles. 

In the past our approach to the determination of the 
unknown lattice spacing of diffraction crystals combined 
two types of crystal lattice spacing measurements: 1) ab- 
solute lattice parameter measurements in which the lat- 
tice parameter of a particular Si crystal is compared to 
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TABLE III: Crystals, crystal orders, and nominal Bragg angles for the various energies. 



Nuclide 


Energy 


Crystals 




A crystal 






B crystal 






(keV) 




order-m 


Brag 


;g angle" 


order-n 


Brag 


;g angle" 


"^Si 


3539 


ILL4.4&6.9 


1 




0.052 


2,-2 




0.105 




3539 


ILL4.4&6.9 


1 




0.052 


3,-3 




0.157 




4934 


ILL4.4&6.9 


1 




0.037 


1 




0.037 




4934 


ILL4.4&6.9 


1 




0.037 


2,-2 




0.075 


33 s 


841 


ILL2.5 


1 




0.220 


1 -1 




0.220 




841 


ILL2.5 


1 




0.220 


3,-3 




0.660 




2380 


ILL2.5 


1 




0.078 


1 




0.078 




2380 


ILL2.5 


1 




0.078 


2,-2 




0.155 




2380 


ILL2.5 


1 




0.078 


-3 




0.233 




5421 


ILL4.4&6.9 


1 




0.034 


1 




0.034 




5421 


ILL4.4&6.9 


1 




0.034 


2,-2 




0.068 




5421 


ILL4.4&6.9 


1 




0.034 


3,-3 




0.102 


36 Qj 


517 


ILL2.5 


2 




0.716 


2 -2 




0.716 




517 


ILL2.5 


2 




0.716 


3,-3 




1.074 




786 


ILL2.5 


1 




0.235 


1,-1 




0.235 




1951 


ILL2.5 


1 




0.095 


2,-2 




0.190 




1951 


ILL2.5 


2 




0.190 


2,-2 




0.190 




6111 


ILL2.5, ILL4.4&6.9 


1 




0.030 


1 




0.030 




6111 


ILL2.5, ILL4.4&6.9 


1 




0.030 


2 




0.061 




6111 


ILL2.5, ILL4.4&6.9 


1 




0.030 


3,-3 




0.091 



the wavelength of an ^^^l2 stabilized laser and 2) lattice 
comparison (relative) measurements in which the small 
lattice spacing difference between known and unknown 
crystal samples was measured. Absolute lattice param- 
eter measurements have been published by researchers 
at the Physikalisch-Technische Bundesanstalt (PTB) in 
1981 8], at the Istituto di Metrologia "G. Colonnetti" 
(IMGC) in 1994 9J, and at the National Measurement 
Institute of Japan (NMIJ) in 1997 0. Lattice com- 
parison measurements have been made at the PTB fTT| 
and the National Institute of Standards and Technology 
(NIST) In Ref. lithe above approach applied to 

the determination of the lattice parameter of the ILL2.5 
crystals is described in detail using the data that was 
available in 1999. Although an early 2004 publication 
contains improved absolute lattice parameter measure- 
ments from IMGC and NMIJ the authors have since 
published an erratum and advised us not to use these re- 
sults 1^. 

One of the authors of the 1998 and 2002 CODATA 
Recommended Values of the Fundamental Physical Con- 
stants has recommended an alternate approach for deter- 
mining lattice parameters values for the crystals used in 
these measurements IT3|. In the adjustment of the fun- 
damental physical constants, absolute and relative lattice 
parameter measurements are used in a consistent way to 
arrive at recommended output values. One of the output 
values is the lattice parameter of an ideal single crystal 
of naturally occurring Si free of impurities and imperfec- 
tions. Since the relative lattice parameter measurements 



connecting the ILL2.5 crystals to the PTB, IMGC, and 
NMIJ absolute lattice parameter crystals in Ref. |6| are 
included in the input data, the value of the lattice pa- 
rameter of the ILL2.5 crystal is an unpublished output 
of the adjustment process. In the 2002 CODATA ad- 
justment only the NMIJ absolute lattice parameter 
value was used as input based on preliminary measure- 
ments that were eventually published in Ref. il4i . Be- 
cause these measurements are now known to be in error, 
the best estimate for the lattice parameter of the ILL2.5 
crystal must be taken from the 1998 CODATA adjust- 
ment [13 and is d(220) ILL2.5 = 1.920155822(57) x IQ-i" 
m at "i? = 22.5 °C in vacuum (relative uncertainty of 
3.0 X 10^^). We prefer to arbitrarily increase the rela- 
tive uncertainty to 5.0 x 10~® to account for the present 
inconsistency of absolute and relative lattice parameter 
results and the variation of the lattice spacing within the 
raw material from which the crystals are manufactured. 

To obtain a value for the lattice parameter of the 
ILL4.4&6.9 crystal, the NIST lattice comparison facil- 
ity was used to measure the lattice spacing difference 
between the ILL2.5 and ILL4.4&6.9 crystals. The di- 
rectly measured relative lattice parameter difference is 
(ILL2.5-ILL4.4&6.9)/ILL2.5 = 4.0 x IQ-^. By combin- 
ing this relative lattice parameter measurement with the 
above absolute value for d(220) ILL2.5 yields the value 
for d(220) ILL4.4&6.9 that is given in Table |IV| The 
reasons for preferring this approach for determining lat- 
tice parameter values over the approach that was used in 
Ref. are all lattice parameter measurements included in 
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TABLE IV: Lattice spacing of ILL2.5 and ILL4.4&6.9 crys- 
tals. The ILL2.5 value is an unpublished output of the CO- 
DATA adjustment process. The NIST lattice comparison fa- 
cility was used to measure the lattice spacing difference be- 
tween the ILL2.5 and ILL4.4&6.9 crystals. 



d(220)" ILL2.5 (m) 


1.920155822(96) x 10'"' 


(ILL2.5 - ILL4.4&6.9)/ILL2.5 


4.0(1.0) X 10"* 


d(220)" ILL4.4&6.9 (m) 


1.920155745(96) x 10""' 


"i? = 22.5 °C, in vacuum 



the 1998 CODATA adjustment are used in a consistent 
way to obtain a value for d(220) ILL2.5 and only one 
direct lattice comparison is needed to obtain a value for 
d(220) ILL4.4&6.9. 

However, it is instructive to compare the d(220) values 
given in Table Hvl to d(220) values obtained with the pro- 
cedure used in Ref.|6t The relative difference between the 
values for d(220) ILL2.5 given in Table Hvl and in Ref.i^ 
is 5.2 x 10"^. The deuteron binding energy and the neu- 
tron mass values given in Ref.^must be corrected for this 
change. In order to use the Ref. approach to determine 
the ILL4.4&6.9 lattice spacing, the ILL 4.4&6.9 crystal 
was compared to the absolute lattice parameter crystals 
from PTE, IMGC, and NMIJ. The relative difference be- 
tween the value for d(220) ILL4.4&6.9 given in Table Hvl 
and the value that is obtained using the Ref. procedure 
is 3.9 X 10~^. These relative differences provide further 
justification for expanding the relative uncertainty of the 
lattice parameter measurements to 5 x 10~®. In addi- 
tion, the lattice parameter difference between the ILL2.5 
and ILL4.4&6.9 crystals can be inferred from the two 
d(220) values obtained with the Ref. 161 approach. The im- 
plied value of (ILL2.5-ILL4.4&6.9)/ILL2.5 = 3.0 x 10"^ 
agrees very well with the directly measured value given 
above and provides further confidence in the lattice pa- 
rameter results that are used to obtain binding energy 
measurements. 



D. Profile recording 

Gamma-ray profiles (intensity vs interferometer 
fringes) were recorded for the crystal configurations given 
in Table ITTll bv scanning the angular position of the sec- 
ond crystal. The profiles consisted of 30 to 45 points with 
counting times from 30 to 180 seconds per point and were 
scanned in both the cw and ccw directions. For each data 
point the interferometer fringe value was reduced to stan- 
dard atmospheric conditions (pressure, temperature, and 
humidity) ^18, 19J. The gamma-ray counts were accu- 
mulated in a Ge detector-multichannel analyzer system. 
The recorded profiles were least squares fit to theoretical 
dynamical diffraction profiles broadened with a Gaussian 
function to account for crystal imperfections, vibrations, 
and thermal and recoil induced motions of the atoms in 




-289.20 -289.15 -289.10 -289.05 364.60 364.65 364.70 364.75 
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FIG. 2; Two representative third order CI profiles along 
with fitted curves: 6111 keV (1,-3) and (1,3) made with the 
thick crystals ILL4.4&6.9. The most important fit parameter 
is the profile centroid. The difference in profile centroids is 
proportional to the Bragg angle while the numeric mean of 
the two centroids gives the offset angle between the second 
crystal diffracting planes and the angle interferometer. 



the source. The use of a single Gaussian function to 
account for deviations of the recorded profiles from the 
theoretical dynamical diffraction profiles has been shown 
to provide reliable peak positions, but is not sufficient to 
obtain nuclear level lifetimes and recoil velocities from 
profile width measurements The adjustable param- 
eters in the fit are the position, intensity, background, 
and Gaussian width contribution. In the fit, the number 
of counts at fringe value i, rii, is weighted by (1/y^)^. 
The detector and fitting procedure are described in more 
detail in Ref. 5. Two representative profiles are shown in 
Fig. 13 The profile positions (in interferometer fringes), 
the most important parameter for the determination of 
energies, are converted into diffraction angles by using 
the second-axis angle calibration that is described in the 
next section. 



IV. ANGLE INTERFEROMETER 
CALIBRATION 

A complete description of a calibration run can be 
found in Ref. "s*. As discussed there the formula connect- 
ing optical fringes /, which are recorded with the rocking 
curves, and the true interferometer angle 9 is 

f = Ksm0 + fo, (2) 

where K is the instrument calibration constant and /o 
is an electronic offset. Both terms must be determined 
experimentally through a calibration procedure. Ideally 
K would be invariable, however experience shows that K 
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depends on temperature, time, humidity and interferom- 
eter laser and alignment. 



Between September 1997 and May 2002, a period 
spanning the measurements described in this paper, the 
GAMS4 spectrometer was calibrated 29 times. Table IVI 
lists each of these calibrations. The experimentally de- 
termined calibration constants (column 5) fit well to a 
linearized equation 

K = Ko + K.g{§ - 26) + Ka{d- 800) +Kh{h- 0.35) 

[-1 ifd< 700, 

+i^lascr X <^ (3) 

[l ifd>700 

where K = K{-d, d, h) is the desired calibration constant, 
'd is the corner cube arm temperature (°C), h is the rel- 
ative humidity, d is the number of days after 8/31/1997, 
and d = 700 corresponds to 8/1/1999 when the interfer- 
ometer laser was replaced. A least squares fit to the data 
in Table IVl gives 

Ko 5133462.12 ±0.59 
Ki) = 7.5 ± 1.3 /°C 
Kd = -0.0059 ±0.0013 
Kh = 41.6 ±3.9 
Xlascr = 3.13 ±0.57 (4) 

where iJavg — 26 °C, h^vg — 0.35, and davg = 800 
(11/9/1999) are average conditions around which the fit 
is made. This fit is plotted in Fig. O The relative stan- 
dard deviation of the residuals (column 7) is 0.33 x 10^^. 
Eqs. (0), 121, and Q along with values for the profile cen- 
troid in fringes, the mean corner cube arm temperature, 
the mean humidity, and the mean wall time can be used 
to extract a profile centroid in radians from a single data 
file. This procedure uses all available calibration data to 
determine one set of coefficients {Kq, K^, K^, Kh, -?Ciasoi) 
that are assumed valid for all of the data presented in 
this paper. In the remainder of this paper this procedure 
is referred to as the Global Calibration. 

The dependence of K on temperature and time is not 
unexpected as the physical dimensions of the interferom- 
eter corner cube arm can vary with temperature and time 
in roughly the amounts seen |^ . The data reveals 
no dependence of K on atmospheric pressure which is 
as expected since the interferometer fringe values are re- 
duced to standard atmospheric conditions and the small 
pressure changes are not likely to alter the dimensions 
of the interferometer. The relative magnitude of the 
dependence on interferometer laser (« 1.22 x 10^^) is 
somewhat larger than the expected variability in the laser 
wavelength and is likely due to interferometer alignment. 
The dependence of K on humidity is unexpected and 
no definitive cause has been established. In an effort to 




12 3 4 



(K,easured-Ko)/Ko (PPM) 

FIG. 3: Global fit to 29 GAMS4 calibrations. Predictions 
from a five parameter global fit [Eq. ^] are plotted ver- 
sus measured calibration values. The standard deviation of 
the relative residuals given in the last column of Table IVI is 
0.33 X 10"^. The dotted lines correspond to ±0.33 x 10"'^ 
about the fit. Across the data sets the temperature ranges 
from 25.5 to 26.8 °C corresponding to a 1.9 x 10~'' relative 
change in K, the relative humidity ranges from 0.21 to 0.58 
corresponding to a 3.0 x 10~^ relative change in K, and 4.6 
years separates the first and last measurements corresponding 
to a 2.0 x 10~® relative change in K. The size of the plotted 
points is proportional to the humidity while darker (lighter) 
points represent lower (higher) temperatures. 



better understand the dependence of the calibration on 
the environment, more frequent calibrations were per- 
formed as the measurements progressed. As is evident 
from Fig. 13 there are significant differences between some 
values of if measured and if fit • Thcse differences prompted 
us to consider alternate calibration procedures. 



B. Other Calibration Procedures 

For the later measurement cycles, multiple calibrations 
were performed: ^^S 1998 - 4, ^^S 1999 - 3 and ^^Si 2000 - 
5. Since these calibrations were dispersed throughout the 
measurement cycle, it is possible to interpolate between 
the measured calibrations using spline fitting to obtain 
calibration constant values for converting profile data in 
fringes to angles. This procedure assumes that each of 
the recorded angle calibrations is valid and the angle cal- 



The Global Calibration Procedure 
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TABLE V: GAMS4 angle calibrations. Each line represents a unique calibration of the spectrometer. Column 1 gives the date 
of the calibration, column 2 the nuclide being measured, column 3 the calibration-average corner cube arm temperature (°C), 
column 4 the calibration-average relative humidity, column 5 the measured instrument calibration constant K, column 6 the 
measured calibration zero, and column 7 the residuals from a five parameter fit to K. The horizontal line indicates the time 
when the interferometer laser was replaced causing a (1.22 ±0.22) x 10~® one time fractional shift in K. "CC" denotes corner 
cube. 





Nuclide 


Avg CC arm 


Avg relative 






Fit 


Date 


measured 


temp (°C) 


humidity 


K 


/o 


rel devxlO*^ 


09/27/1997 


36 


25.595 


0.495 


5133471.3 


588.3 


-0.92 


03/08/1998 




26.272 


0.317 


5133463.2 


-747.7 


0.01 


03/26/1998 




26.266 


0.209 


5133460.1 


-721.1 


-0.28 


03/30/1998 




26.256 


0.308 


5133462.7 


-721.2 


-0.01 


09/26/1998 


33s 


26.716 


0.481 


5133473.7 


-547.5 


-0.30 


10/21/1998 


33 s 


26.607 


0.323 


5133461.8 


-547.7 


0.55 


11/05/1998 


33 s 


26.308 


0.314 


5133462.6 


-203.1 


-0.11 


11/17/1998 


33s 


26.199 


0.265 


5133458.1 


-202.9 


0.18 


04/13/1999 


33s 


26.041 


0.323 


5133458.1 


-186.9 


0.25 


04/30/1999 


33 s 


26.194 


0.536 


5133466.0 


-186.9 


0.64 


05/07/1999 


33s 


26.173 


0.583 


5133471.2 


-68.2 


-0.02 


09/17/1999 




26.492 


0.464 


5133474.4 


-268.2 


-0.09 


09/19/1999 




26.486 


0.434 


5133473.6 


-268.1 


-0.18 


09/25/1999 




26.648 


0.503 


5133474.1 


-255.5 


0.52 


10/05/1999 




26.354 


0.334 


5133468.8 


-254.5 


-0.26 


10/15/1999 




26.348 


0.449 


5133470.3 


-249.5 


0.37 


10/18/1999 




26.477 


0.369 


5133469.3 


-249.4 


0.08 


12/03/1999 




25.493 


0.257 


5133454.5 


-249.5 


0.59 


07/06/2000 




26.545 


0.399 


5133471.3 


-183.5 


-0.26 


07/18/2000 




26.474 


0.390 


5133468.7 


-183.4 


0.05 


10/15/2000 


^'^Si 


26.779 


0.412 


5133471.0 


-253.2 


0.13 


10/22/2000 


^'^Si 


26.490 


0.374 


5133468.5 


-192.4 


-0.12 


11/05/2000 


^'^Si 


26.313 


0.314 


5133465.4 


-240.8 


-0.28 


11/14/2000 


^'^Si 


26.409 


0.320 


5133464.8 


-262.2 


0.00 


11/28/2000 


29Si 


26.302 


0.312 


5133462.9 


-276.1 


0.14 


10/28/2001 




26.725 


0.374 


5133469.1 


-211.8 


-0.32 


11/06/2001 




26.619 


0.316 


5133465.0 


-250.0 


-0.16 


11/12/2001 




26.347 


0.236 


5133459.4 


-250.0 


-0.13 


05/19/2002 




26.551 


0.411 


5133466.7 


-263.4 


-0.05 



ibration varies smoothly with time during the measure- 
ment cycle. Although the individual calibration values 
within a measurement cycle are dependent on tempera- 
ture, humidity, and day number, the spline fitting uses 
only the individual calibration values and assumes no ex- 
plicit dependence on temperature, humidity or day num- 
ber. We have called this calibration procedure the Spline 
Calibration. It uses only the calibration information ob- 
tained during a given measurement cycle to analyze the 
wavelength data recorded in that cycle. The ■^^S and ^^Si 
data were analyzed using the spline calibration and the 
relative difference between wavelengths obtained with the 
global and spline calibration procedures is « 0.2 x 10~^. 
Since only one calibration was recorded during the ^^Cl 
1997 measurement cycle, the spline calibration could not 



be applied to this data. This led us to a third calibration 
procedure called the Local-Global Calibration. 

In the local-global calibration, the calibrations that 
were performed in a particular measurement cycle are 
given a more significant role in the determination of 
K{{},d,h) for that measurement cycle. First, K^, Kd, 
Kh, and if laser are taken as given in Eq. |0J since these 
coefficients are not expected to change with time and a 
global fit provides the best estimate of these coefficients. 
Next the subset of calibrations performed in each mea- 
surement cycle are fit to the equation 

K= Xo + 7.5(?9- 26)- 0.0059(d- 800) 

+41.6(/^ -0.35) + 3.13 x/-^ ^ < 700, 

1 if d > 700 
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to obtain independent values of Kq for each measure- 
ment cycle. All of the ^^S and ^^Si data were analyzed 
using the local-global calibration procedure and the rela- 
tive difference between wavelengths obtained with global 
and local-global calibration procedures is « 0.3 x 10^^. 
The "^^Cl data were also analyzed using the local-global 
calibration procedure and show significantly larger rela- 
tive differences (~ 0.9 x 10^^) between wavelengths ob- 
tained with the global and local-global calibration proce- 
dures. These large differences result from recording only 
one calibration during the '^^Cl measurement cycle and 
from the fact that this particular calibration is the most 
discrepant point in the global fit of Fig. |3| 

Consideration of these alternate calibration procedures 
did not provide sufficient evidence to make a "best" cal- 
ibration procedure choice. All of the wavelength and en- 
ergy values reported in this paper were obtained using 
the global calibration procedure. Our reasons for choos- 
ing this procedure are 1) all of the available data can 
be analyzed using one procedure and 2) this procedure 
makes the maximum use of the available calibration data. 



C. Calibration Uncertainty 

Although consideration of three calibration procedures 
did not lead to a clear "best" procedure, this exercise 
does provide an estimate of the calibration uncertainty. 
The variation of wavelengths obtained with the global, 
spline and local-global calibration procedures suggests a 
relative calibration uncertainty of 0.2-0.3 x 10~®. A sec- 
ond estimate of the calibration uncertainty is available 
from the fit used in the global calibration procedure. The 
standard deviation of the relative residuals given in Ta- 
ble El a-nd shown in Fig. O is 0.33 x 10~^ and provides 
a measure of the quality of the fit and the uncertainty 
of this calibration procedure. A third estimate of the 
calibration uncertainty is the variation of the wavelength 
values obtained in different measurement cycles. In this 
approach lower energy (larger angles) intense transitions 
must be used because higher energy (smaller angles) weak 
transitions have a statistical uncertainty that masks the 
calibration uncertainty. The 841 keV line in '^^S was mea- 
sured in 1998 and 1999 and shows a relative excess vari- 
ation of 0.38 X 10~^ if the global calibration procedure is 
used and 0.16 x 10~^ if the spline calibration is used. In 
addition, this instrumental set up was used to measure 
an intense line that does not contribute to the binding 
energy determinations presented here, namely the 816 
keV line in ^^*Er. This line was measured in two dif- 
ferent measurement cycles, Oct-Nov 2000 and Nov 2003 
and shows a relative excess variation of 0.45 x 10~^ and 
0.3 X 10^^ for the global and spline calibrations, respec- 
tively. Although it is difficult to obtain a rigorous cali- 
bration uncertainty from these three estimates, we choose 
to use the above values to arrive at a slightly conserva- 
tive relative calibration uncertainty of 0.4 x 10~^. This 
calibration uncertainty will be combined with the statis- 



tical and other systematic uncertainties to obtain final 
wavelength uncertainties. 

V. WAVELENGTH MEASUREMENTS 

Wavelengths are determined by combining a sequence 
of profile angle measurements. First the profile centroids 
and uncertainties in interferometer fringes are converted 
into angles via Eq. Q . The calibration constant K that 
appears in this equation has been appropriately corrected 
for the corner cube arm temperature, the relative hu- 
midity, and time. Each of these angles is next corrected 
for vertical divergence as discussed in Ref. I^J Typically 
these corrections are between 2 and 4 x 10"'' in relative 
size. For a given energy and configuration the profiles are 
recorded with the first crystal in a fixed position and the 
second crystal sequentially in a more and less dispersive 
position (see Fig. In addition, the data recording se- 
quence includes both cw and ccw rotation of the second 
crystal. To more concretely illustrate how wavelengths 
are determined, we use the ^^S 841 keV (1,-3), (1,3) 
measurement as an example. A group of four profiles 
recorded in the sequence (1,-3 cw), (1,3 cw), (1,3 ccw), 
(1,-3 ccw) is used to determine a wavelength value. The 
four angles associated with the four profiles are fit with 
the equation 

0(n,i,^)=arcsin(^^^^)+0o(i), (6) 

where n is the diffraction order, t is the time, is the 
crystal temperature, Amoas is the sought after wavelength, 
d^t}) is the lattice spacing at crystal temperature i5 and 
Ooit) represents the potentially time dependent angular 
offset between the second crystal diffracting planes and 
the angle interferometer. The symbol Amoas is introduced 
to indicate that wavelengths determined using this equa- 
tion are laboratory measured wavelengths (not corrected 
for recoil); d{'&) is given by 

dm = ^22.5 °c, atm(l + 2.56 X 10-^(i9 - 22.5)) , (7) 

where ^22. 5 °c, atm is the lattice spacing at 22.5 °C and 
atmospheric pressure. The lattice parameter measure- 
ments given in Table IIVI are specified for vacuum. To 
obtain the value at the pressure present in the reactor 
hall {p w 0.987 atmospheres) it is necessary to use the 
following transformation 

^22.5 "dp) = C?22.5 °C, vac(l " ^P) , (8) 

where p is the pressure. In this equation e — 0.3452 x 
lO-'^/atmosphere [HEi- 

In the fit to Eq. (O, the three parameters are the wave- 
length, a constant angular offset, and a linear temporal 
offset term (i.e., 6o{t) = a+bt). Each of the four angles di 
is weighted by (l/cej^ where ag^ is the profile centroid 
uncertainty. This sequence of profiles is repeated multi- 
ple times in at least two different sets of orders (here 29 



9 



instances of (1,-3), (1,3); 21 instances of (1,-1), (1,1); 
and 1 instance of (1,-3), (1,-1), (1,1), (1,3)). An uncer- 
tainty equal to the standard deviation of the wavelength 
determinations in a unique set of orders is assigned to 
each wavelength determination derived from that set of 
orders. The average wavelength is the weighted mean of 
all the individual determinations. We find no evidence 
for an order dependent effect in the wavelength data. 

Table I VII gives the nuclide, the nominal energy, the 
crystals, the total number of Bragg angle measurements 
and mean wavelength along with the statistical uncer- 
tainties in parentheses. The sulfur transitions appear 
twice because they were measured in 1998 and 1999. 

In Table I VI II final measured wavelength values and 
uncertainties are reported (column 3). Four additional 
sources of uncertainty have been added in quadrature to 
the statistical uncertainties given in Table IVTI These are 
the calibration uncertainty discussed in Sec. II VI and un- 
certainties associated with the crystal temperature, the 
vertical divergence of the gamma-ray beam, and the mea- 
sured lattice spacing. A relative calibration uncertainty 
of0.4xl0~^is applied to each energy listed in Table IVll 
To obtain final sulfur wavelengths and uncertainties it 
is necessary to add the calibration and statistical uncer- 
tainties in quadrature for each sulfur transition before 
statistically combining the 1998 and 1999 values. This 
has the effect of reducing the calibration uncertainty for 
the sulfur lines because they were measured twice. In the 
case of binding energies where one sums transition ener- 
gies, the calibration uncertainty is applied to the sum 
rather than to the constituent transitions; and again, 
two binding energies are combined to arrive at a final 
sulfur binding energy. The other three uncertainties are 
applied to the final transition energies or to the binding 
energies. An uncertainty in the crystal temperature mea- 
surement of 0.05 °C contributes a relative uncertainty of 
0.1 X 10~^. The vertical divergence uncertainty accounts 
for the possible misalignment of the gamma-ray beam 
with respect to the plane of dispersion of the spectrom- 
eter. A misalignment of 5 mm over a distance of 15 m 
leads to a relative wavelength uncertainty of 0.06 x 10~^. 
This uncertainty is approximately 20% of the size of the 
correction. From Sec. IIII CI and Table IIVI the relative 
crystal lattice spacing uncertainty is 0.05 x 10~^. 



VI. BINDING ENERGY MEASUREMENTS 

In column 4 of Table IVIII the corresponding energy 
equivalent values i?moas are given, where the conversion 
factor Im-i = 1.23984191(11) x lO"*' eV was used ^2]. 

The measured wavelength values must be corrected 
for recoil to obtain wavelength values Atrans whose cor- 
responding energies can be summed to obtain binding 
energies Sn- To excellent approximation we account for 



the recoil by the transformation 



Atrans Anicas lA'l Cp" y Anicas / 

where M is the mass of the decaying nucleus, h is the 
Planck constant, and c is the speed of light. The sec- 
ond term in Eq. (j5J accounts for the loss of energy im- 
parted to the recoiling nucleus. If the decay occurs in 
flight, there will be a first order Doppler effect. It causes 
no shift in the central value as long as the motion is 
isotropic as is expected to be the case here. As discussed 
in Ref. l2fiL two additional terms appear in Eq. © in the 
case of decays from the capture state. First, the kinetic 
energy of the incident neutron (« 0.057 eV) must be sub- 
tracted from the measured gamma-ray energy. For the 
three capture gamma-rays measured here (3539 keV in 
29Si, 5421 keV in ^Ss, and 6111 keV in ^sq), this effect 
is less than 0.02 x 10~^. Second, there is a Doppler term 
if the incident neutron comes from a particular direc- 
tion. The relative peak-to-peak amplitude of this term is 
K, 0.6 x 10^® for the nuclei being discussed here. The term 
vanishes if the incident neutron direction is isotropic as 
is expected to be the case here. The relative uncertainty 
of the recoil correction is estimated to be no greater than 
0.01 X 10~^. As such it is negligible given the current ac- 
curacy of Anioas. Values for Atrans and the corresponding 
energy equivalent i?trans are given in Table IVm column 
5 and 6. 

To obtain a wavelength equivalent to the binding en- 
ergy, Abo, we sum the reciprocals of the Atrans values for 
the transitions comprising the binding energy cascade 
and then convert Abe into atomic mass units using the 
conversion factor: Im^^ = 1.3310250506(89) x lO'^^ u 
0. Likewise, the binding energy can be expressed in eV 
by using the m~^ to eV conversion factor given above. 
Table IVIIll contains values for the three binding energies 
of ^^Si, '^^S, and '^^Cl in meters, atomic mass units and 
electron volts. 



VII. DISCUSSION 

In this section we discuss the consistency of gamma- 
ray based and atomic mass based binding energy mea- 
surements and compare the gamma-ray measurements 
reported here with other high precision gamma-ray mea- 
surements. As discussed in Sec. 13 precision atomic mass 
measurements can be used to determine binding energies. 
By usin g at omic mass values from the Atomic Mass Data 
Center [271 l2a | and the fundamental constants 2] along 
with Eq. (^) , binding energies primarily based on atomic 
mass measurements can be determined. However, since 
the determination of the neutron mass includes a gamma- 
ray measurement, binding energy determinations based 
on Eq. ^ are a combination of atomic mass measure- 
ments and gamma-ray measurements. This difficulty can 
be circumvented by expressing m(n) in terms of m(H), 
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TABLE VI: Values of the measured wavelengths for the transitions included in the binding energy determinations. The 
uncertainties are statistical only. 



Nuclide 


Energy 


Crystals 


Number of 


■^mcas X 10 


relative 




(keV) 




Bragg angle 
measurements 


(m) 


uncertainty 
xlO*^ 


^^Si 


3539 


ILL4.4&6.9 


42 


0.350340126(44) 


0.13 




4934 


ILL4.4&6.9 


147 


0.25128808(18) 


0.70 


33s 


841" 


ILL2.5 


26 


1.47429306(14) 


0.09 




2380" 


ILL2.5 


18 


0.52103852(30) 


0.58 




5421" 


ILL4.4&6.9 


42 


0.228730970(63) 


0.27 




841' 


ILL2.5 


25 


1.47429225(12) 


0.08 




2380' 


ILL2.5 


31 


0.52103905(23) 


0.43 




5421' 


ILL4.4&6.9 


30 


0.22873089(13) 


0.56 


36 Q 


517 


ILL2.5 


15 


2.39782393(10) 


0.04 




786 


ILL2.5 


2 


1.57681233(83) 


0.53 




1951 


ILL2.5 


12 


0.63544928(10) 


0.16 




6111 


ILL2.5, ILL4.4&6.9 


17 


0.20288757(10) 


0.50 



"Sept-Oct 1998 
'April-May 1999 



TABLE VIL Measured and recoil corrected (transition) wavelengths and energies. Measured values from Ref. 0for ^H, which 
have been adjusted for the adjusted value of d(220) ILL2.5, are included for convenience. 



Nuclide 


Energy 
(keV) 


-^mcas X 10 

(m) 


-£'mcas 

(eV) 


Atrans X 10"*^^ 

(m) 


-fi^trans 

(eV) 


wavelength relative 
uncertainty 
xlO'* 


^^Si 


3539 
4934 


0.35034013(15) 
0.25128808(21) 


3538966.3(1.6) 
4933946.3(4.0) 


0.35031716(15) 
0.25126511(20) 


3539198.3(1.6) 
4934397.4(4.0) 


0.44 
0.82 


33s 


841 
2380 
5421 


1.47429265(47) 
0.52103883(24) 
0.228730944(95) 


840974.08(28) 
2379557.6(1.1) 
5420525.5(2.3) 


1.47427246(47) 
0.52101864(24) 
0.228710758(95) 


840985.60(28) 
2379649.8(1.1) 
5421003.9(2.3) 


0.32 
0.47 
0.42 


36 Q 


517 
786 
1951 
6111 


2.3978239(10) 
1.5768123(11) 
0.63544928(29) 
0.20288757(13) 


517069.62(22) 
786296.43(53) 
1951126.47(89) 
6110980.2(4.0) 


2.3978054(10) 
1.5767938(11) 
0.63543077(29) 
0.20286906(13) 


517073.61(22) 
786305.66(53) 
1951183.30(89) 
6111537.6(4.0) 


0.42 
0.67 
0.45 
0.66 




2223 


0.557671328(99) 


2223248.44(44) 


0.557341007(99) 


2224566.10(44) 


0.18 



TABLE VIII: Measured binding energies in meters (m), atomic mass units (u) and electron volts (eV). Abe = ^ i - is the 

wavelength of a photon whose energy is equal to the binding energy. Measured values from Ref. 6 for ^H, which have been 
adjusted for the adjusted value of d(220) ILL2.5, are included for convenience. 



Nuclide 


Abo X 10^2 
(m) 


Sn X 10^ 

(u) 


(eV) 


wavelength relative 
uncertainty 
xlO** 


29Si 


0.146318275(86) 


9.0967793(53) 


8473595.7(5.0) 


0.59 


33s 


0.143472991(54) 


9.2771820(35) 


8641639.8(3.3) 


0.38 


36 Qi 


0.144507180(80) 


9.2107883(51) 


8579794.5(4.8) 


0.55 




0.557341007(98) 


2.38816996(42) 


2224566.10(44) 


0.18 
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m(^H) and S'„(^H). By making this substitution and re- 
arranging terms, Eq.^ becomes 

m(^X) - m(^+iX) + m(2H) - m(H) = 

S,,{^+'X) ~ SnCR) . (10) 

Since the left and right sides of this equation involve only 
atomic mass and gamma-ray measurements respectively, 
this equation is a valid test of the consistency of high- 
precision atomic mass and gamma-ray measurements. 

Recently, new highly accurate values for the left hand 
side of this equation have been reported for ^"""^X equal 
to ^^Si and '^^S In these measurements the cyclotron 
frequencies of two different ions simultaneously confined 
in a Penning trap were directly compared. The measured 
quantities are the mass ratios to['^^S'^]/to['^^SH'^] and 
m[^^Si~'~]/m[^®SiH^] from which, along with the quan- 
tity to(^H) — 2to(H), the mass differences on the left hand 
side of Eq.^lare derived. In Table Hxl column 2 the left 
hand side mass differences for Si and S are given. The 
relative uncertainties for the Si and S values in column 2 
are 7.0 x 10~® and 7.3 x 10~* respectively. 

Values for the right hand side of this equation for '^^^X 
equal to ^^Si and ^^S follow directly from the binding en- 
ergies given in Table IVlTTl and are given in Table Hxl col- 
umn 3. The relative uncertainties for the Si and S values 
in column 3 are 8.0 x 10^*" and 5.1 x 10~^ respectively. 

In Fig. 2] the values given in Table IIXI are plotted to 
show the consistency of the atomic mass and gamma- 
ray measurements. The left and bottom axes are used 
for ^^Si, while the top and right axes are used for ^^S. 
The scales of the axes have been chosen so that a diago- 
nal line through the plot represents exact consistency be- 
tween atomic mass and gamma-ray measurements. The 
plot shows that the quality of the consistency test is lim- 
ited by the uncertainty of the gamma-ray measurements. 
For ^^Si the two measurements are slightly inconsistent 
(1.2 cr), while for ^^S the two measurements agree within 
the uncertainty (0.4 a). 

In Table llXl column 4 the fractional difference between 
the atomic mass and gamma-ray measurements is given 
along with the weighted average of the Si and S fractional 
differences. These measurements confirm the consistency 
of atomic mass and gamma-ray measurements with a rel- 
ative uncertainty of 4.3 x 10~^. 

The most accurate Si and S gamma-ray energies were 
measured using Ge(Li) solid state spectrometers. These 
spectrometers derive an energy scale from gamma-ray 
standard energies. For high precision comparisons, the 
published values of the gamma-ray energies need to be 



adjusted to account for changes in the energy standards. 
Because a number of energy standards covering a wide 
energy range are used, the shift in the energy scale for a 
particular energy is difficult to estimate. Although the 
procedure that has been used is not rigorous, it is likely 
sufficient given the accuracy of the published energies. 
New values for the standards were taken from Ref. |33 
or have been determined using atomic mass values from 
Refs.HIH andlH For ^^Si the values of the 3539 keV 
and the 4945 keV energies in Ref. |^ have been corrected 
by — 7 eV to account for the change in the 2223 keV and 
the 4945 keV standards produced in the "'^H(n,7) and 
^^C(n,7) reactions, respectively. For ^^S the values for 
the 841 keV, 2380 keV, and 5421 keV energies in Ref. HI 
have been corrected by -1-8 eV, —107 eV, and —107 eV re- 
spectively. These corrections account for changes in the 
412 keV standard produced in the decay of ^^*Au and 
changes in the 2223 keV, the 4945 keV, and the 10829 
keV standards produced in the ^H(n,7), the ^^C{n,'-f), 
and the ^^N(n, 7) reactions, respectively. For the low en- 
ergy CI lines considerably more precise data exists. In 
1985 energy values for some CI lines were measured us- 
ing the GAMS4 facility in its early stage of development 
|33| . These published values also need to be corrected 
for changes in the fundamental constants and known 
errors in the lattice spacing of the crystals ^3^. These 
corrections to the CI lines can be made with much more 
certainty than the corrections to the ^^Si and '^^S lines. 

In Tabled we compare the -Emcas values in this report 
(column 2) with the corrected -Emcas values from other 
references (column 3). For the ^^Si and ^^S gamma-rays 
the new measurements agree with the corrected older 
measurements within the uncertainty except for the 5421 
keV line which differs by 1.3 times the combined uncer- 
tainty. Because of the large uncertainty of the older Si 
and S measurements, this comparison does not provide a 
very stringent test of our new measurements. However, 
the consistency of the CI measurements over more than 
15 years during which the spectrometer, the crystals, and 
measurement procedures were significantly changed lends 
a large measure of confidence to the gamma-ray measure- 
ments. 
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TABLE IX: Consistency of high precision atomic mass and gamma-ray measurements. 





m('*X) - m('*+^X) + m(^H) - m(H) 


S^f-^+^X) - Sr.C^R) 


rGlcitive difference 




(u) 


(u) 


(col3 - col2)/col3 X 10^ 




0.00670861569(47) 


0.00670860929(536) 


-9.54(8.02) 


33s 


0.00688901053(50) 


0.00688901206(351) 


2.22(5.15) 


weighted average 








relative difference 






-1.21(4.33) 



TABLE X: Comparison of measured gamma-ray energies. 



Nuchde 


This report 


Other references" 




F, 

J-^mcas 


-£^mcas 




(eV) 


(eV) 


2^Si 


3538966.3(1.6) 


3538973(40) 




4933946.3(4.0) 


4933973(30) 


33s 


840974.08(28) 


840982(14) 




2379557.6(1.1) 
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FIG. 4: Mass differences determined in a Penning trap versus 
those determined from wavelength determinations. The left 
(right) and bottom (top) axes correspond to Si (S). The scales 
of the axes have been chosen so that a diagonal line through 
the plot represents exact consistency between Penning trap 
and gamma-ray measurements. The plot shows that the qual- 
ity of the consistency test is limited by the uncertainty of the 
gamma-ray measurements. 
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